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Abstract 

The quantum stochastic Schrodinger equation or Hudson-Parthasareathy (HP) equa¬ 
tion is a powerful tool to construct unitary dilations of quantum dynamical semigroups 
and to develop the theory of measurements in continuous time via the construction of 
output fields. An important feature of such an equation is that it allows to treat not only 
absorption and emission of quanta, but also scattering processes, which however had very 
few applications in physical modelling. Moreover, recent developments have shown that 
also some non-Markovian dynamics can be generated by suitable choices of the state of 
the quantum noises involved in the HP-equation. This paper is devoted to an application 
involving these two features, non-Markovianity and scattering process. We consider a 
micro-mirror mounted on a vibrating structure and reflecting a laser beam, a process giv¬ 
ing rise to a radiation-pressure force on the mirror. We show that this process needs the 
scattering part of the HP-equation to be described. On the other side, non-Markovianity is 
introduced by the dissipation due to the interaction with some thermal environment which 
we represent by a phonon field, with a nearly arbitrary excitation spectrum, and by the 
introduction of phase noise in the laser beam. Finally, we study the full power spectrum 
of the reflected light and we show how the laser beam can be used as a temperature probe. 

Keywords: Quantum optomechanics, quantum stochastic differential equations, radiation 
pressure interaction, quantum Langevin equations, heterodyne detection. 

1 Introduction 

Quantum optomechanical systems represent an active field of research, very important both 
from the theoretical and experimental point of views, with applications in quantum optics 
and quantum information A great interest is due to the possibility of seeing quantum 

effects in a macroscopic mechanical resonator, say a mirror mounted on a vibrating structure 
and coupled to optical elements by radiation pressure. Typically, the theoretical description of 
such a kind of systems is based on the quantum Langevin equations 171181 . a flexible approach 
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allowing also for the introduction of non Markovian effects. Recently, some experimental 
evidence of non Markovian effects in an optomechanical system has been reported ||9l- 

To get mathematically consistent quantum Langevin equations one has to use the quantum 
stochastic calculus and the quantum stochastic Schrodinger equation, or Hudson-Parthasarathy 
equation (HP-equation) nniin]. In this mathematical context the quantum Langevin equa¬ 
tions appear under the name of Evans-Hudson flows or quantum flows GIHSl. In m a 
description of a dissipative mechanical oscillator has been obtained in terms of quantum 
stochastic differential equations; this description is fully consistent and valid at any temperat¬ 
ure and it respects some symmetry requirements and physical constraints such as a weak form 
of equipartition at equilibrium and the translation invariance of the dissipative part of the dy¬ 
namics. Non Markovian effects have been introduced by a suitable choice of the state of the 
quantum noises appearing in the HP-equation. In Sectionl^the quantum stochastic model of 
a dissipative mechanical oscillator is presented. An equation of Hudson-Parthasarathy type 
gives the unitary dynamics of the oscillator interacting with a Bose held (here representing 
the phonon held). The evolution equations for the system operators in the Heisenberg picture 
are the quantum Langevin equations and a suitable choice of the state of the held (based on a 
held analog of the P-representation in the case of discrete modes) allows for the introduction 
of thermal, non Markovian effects. 

Usual quantum Langevin equations allow to describe absorption and emission of energy 
quanta by the main system, not scattering processes. As a matter of fact it seems that the 
existing literature takes into account the mirror/light radiation pressure interaction only if 
mediated by cavity modes; indeed the subject is often called cavity optmechanics flUlgl. This 
is due to the fact that some interesting physical phenomena as laser cooling appear when 
a cavity mode is involved ll3l-l6l fT4l . but also to the fact that a way to describe at a quantum 
level the direct scattering of laser light by a vibrating mirror is lacking. In this respect, another 
advantage of the HP-equation is that it allows also for the description of scattering processes 
lUSlIIhl. The content of Section [3 will be to introduce in the HP-equation, describing the 
mechanical oscillator, the radiation pressure interaction due to a laser directly illuminating the 
mirror. Finally, in Section|4]we show how to describe the heterodyne detection of the reflected 
light and we study the properties of the resulting power spectrum; this last step involves also 
the theory of measurements in continuous time EHSa. In particular we show that, for a 
weak probe laser, the model produces explicit expressions for the spectrum due to elastic 
scattering of the photons and for the side-bands due to Stokes and anti-Stokes scattering. 

We end this section by introducing the HP-equation. We give a short, heuristic present¬ 
ation; a mathematically rigourous formulation of quantum stochastic calculus, HP-equations 
and related notions can be found in 01 Illl2lfr9ll21ll2^ . Firstly, we introduce the formal helds 
bk{t), t G IS., k = 1... d, satisfying the canonical commutation rules (CCRs) 


bfls),bl{t) 




[bi{s),bk{t)] = 0. 


( 1 ) 


In this paper we consider only the representation of the CCRs O on the Fock space, the 
one characterized by the existence of the vacuum state. For quantum stochastic calculus 
involving non-Fockrepresentations see, for instance, ll24l . Let us introduce the Hilbert space 
= L^(R; C'^) (the one-particle space) and its symmetrized powers 
(the n-particle space). We denote by F = r(L^(]R.; C^)) the symmetric Fock space over 
L^(R; C^), i.e. F = C © J2n=i and by e(/), / G L^{R- C^), the coherent 

vectors, whose components in the 0,1,..., n,... particle spaces are 


e(/) := e-^ (l, /, ■ (2) 
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Note that e(0) represents the vacuum state and that 


(e(5)|e(/)) = exp ||/f - i \\gf + {g\f)^ . 


Let {zfe, fc > 1} be the canonical basis in and for any / € L^(R; C'^) let us set fk{t) '■= 
{zk\ f(t))cd. Then we have bk(t) e{f) = fk{t)e{f). By formally writing 


Bk{t) = 



6fe(s)ds, 



(3) 


we get the annihilation and creation processes, families of mutually adjoint operators, whose 
actions on the coherent vectors are given by 


Bk{t) 


3(/) = / fkis) 
Jo 


dse(/). 


{eig)\Bl(t)eif)) = 


Pfe(s)ds(e(p)|e(/)). 


The overline denotes the complex conjugation. It is a property of the Fock spaces the fact 
that the action on the coherent vectors uniquely determines a densely dehned linear operator. 
In terms of the integrated processes the CCRs ([T]l become 

= 4imm{f,s}, [Bl{t), b}( s)] = [Bkit), Bi{s)] = 0. (4) 


We introduce also the gauge processes 


Akiit)=f bl{s)bi{s)ds, {e{g)\Akiit)e{f)) = f gkis) fi{s)ds {e{g)\e{f)). (5) 

Jo Jo 


The operator Akk (t) turns out to be a number operator and it counts the quanta present in the 
held k in the time interval (0,f). Quantum stochastic calculus is an ltd type calculus with 
respect to the integrators dt, dBk{t), dBl{t), dAki{t) satisfying the ltd product rules 


dBk{t)dBj{t) = Skidt, 

dAki(t)dBl{t) = 6iidBl(t), 


dBi{t)dAki{t) = 6ikdBi{t), 
dAki{t)dAij(t) = 6iidAkj{t); 


(6) 


all the other possible products vanish. We shall need also the generalized Weyl operators 
W{g', V), where g G C^) and 14 is a unitary operator on C^); these are unitary 

operators dehned by 

W(g;V)e(f)=exp{iIm(Vflg)}e(Vf + g), Vf G L^(R;C‘^). (7) 


From the dehnition one obtains the composition law 

W(h; V) yV{g; U) = exp {—ilm {h\Vg)} yV{h + Vg; VU). 
In the case V = 1, it is possible to show that 

p + OO 


( 8 ) 


>V(p;l) =exp|^(^y 5 fc(f)dBj;(f)-h.c.^j, 


where h.c. means hermitian conjugate, and from O one sees that yV(g; 1) is the held analog 
of what is called a displacement operator in quantum optics ll22l . 
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Let us introduce now a quantum system with separable Hilbert space J-C; let H, Rk, Ski 
be system operators with H self-adjoint and the operator matrix (Ski) debning a unitary 
operator S' on J{ (g) C^. Then, we consider the system/held evolution equation given by the 
HP-equation 




k kl 



(9) 


with the initial condition U (0) = 1. When H and Rk are bounded operators there is a 
unique solution, which is unitary and strongly continuous in t mi]. When these operators 
are unbounded some restrictions are needed in order to control the domains; then, existence, 
uniqueness, unitarity can be proved insiisi. The solution U (t) gives the evolution in the 
interaction picture with respect to the free evolution of the held, which is modelled by the so 
called left time shift 0(f) in the Fock space; indeed, U(t) = Q{t)U{t), t > 0, and U(t) — 
[/(—f)^0(—f)^, f < 0, defines a strongly continuous unitary group, whose Hamiltonian has 
been characterized in ESi . 

If we now consider a generic system operator X, its evolution in the Heisenberg descrip¬ 
tion is given by X{t) = U{t)^XU{t). By differentiating this product according to the rules 
of quantum stochastic calculus, summarized by •HJ, and taking into account that U(t) is a 
unitary operator, we get the quantum Langevin equations 



k 


+ Y, Sik{t)Hx{t), Riit)]dBl it) - Yit)^]Sik{t)dBk it) 


kl 


kl 



If po is a generic statistical operator for the system and a a held state, we can consider the 
reduced state of the system p(f) = Trr {C/(f)po ® CTC/(f)i}. When a is the vacuum state or, 
more generally, a coherent vector, then the reduced system state pit) satishes a Markovian 
master equation GIlIIll with a Lindblad type generator Il26l . If a more general state is taken 
for cr, non-Markov effects enter into play and a simple closed evolution equation for the 
reduced dynamics could even not exist 018lll9ll27l . 

Also the helds in the Heisenberg picture can be introduced ll28l : these are the output fields 


Br\t) =Uit)^Bkit)Uit), ^(0 = Uit)^Blit)Uit), 

Airit) = Uit)^Akiit)Uit). 


( 11 ) 


The outputs helds represent the helds after the interaction with the system, while Bkit), 


Blit), Akiit) are the helds before the interaction and, so, they are called input fields. By 


differentiating the products dehning the output helds and using (|3 and (|6]l, we get the in¬ 
put/output relations im 


dBl'^\t) = Y Skiit)dBkit) + Rkit)dt, 


( 12 ) 
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dA°f (t) = ^ SUt)^Riit)dBlit) 

ij i 

Y, +Rk{t)^ Su{t)dB,{t) + Rk{t)^ Ri{t)dt. (13) 

i 

By the properties of U{t) we get U(T)^Bk{t)U(T) = U(t)^Bk{t)U{t), \/T > t, and similar 
equations for the other fields. This implies that the output fields satisfy the same CCRs as 
the input fields. Self-adjoint combinations of the output fields commuting for different times 
represent field observables which can be measured with continuity in time and this is the key 
ingredient for a quantum theory of measurements in continuous time iniiisiEi]. 


2 Langevin equations for a mechanical oscillator in a 
thermal bath 


In this section we present the description of a quantum dissipative mechanical oscillator ob¬ 
tained in ifTrl Sects. 2, 3]. The Hilbert space of the system is TC = L^(R) and q and p denote 
the usual position and momentum operator in dimensionless units, satisfying the commuta¬ 
tion relations [q,p] — i. We denote by > 0 the bare frequency of the mechanical oscil¬ 
lator and by 7i„ > 0 its damping rate; we consider only the underdamped case: > 7m/2. 

Then, we introduce the damped frequency Wm and the phase factor r by 




2 H 

n 


(14) 


We define now the mode operator 

am = \ {q + irp) = ^ (flraq +‘^p + iuJmp) , (15) 

satisfying the commutation rules [am, Q^m] = 1- The inverse transformation turns out to be 

{t am + rai,) , P = ^ (“m “ am) • ( 16 ) 


9 = 

We need also the self-adjoint operator 
hvt 


iTm = 


{p^+q^)+ ^^{q,p} = + ^ ■ (17) 


We introduce now the HP-equation for a mechanical oscillator in a thermal bath by taking 
in ® a single field Bi{t) = i?th(i) and H = H^, i?i = ^ = 1. By (fTOl i the 

quantum Langevin equations for a^n, q,P turn out to be 


dain(f) = - (iWm + ain(f)df - y/^dBth{,t), 


( 18 ) 


dq{t) = n„,p{t)dt -I- dC',(f), 

dp{t) = -(n^qit) + 'Y^p{t))dt + dCp{t), 


(19a) 

(19b) 
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in which we have introduced the Hermitian quantum noises 


c,{t) = -\/H^. 

By (|4|i the new noises obey the commutation rules 


(20) 


[Cq{t),Cp{s)] = i7mmin{t,s}, [Cq{t),Cq{s)] = [Cp{t),Cp{s)\ = 0. (21) 

Obviously, from (fT5l l. (fTSl l we have that the equation (fTSl) for Om is equivalent to the system 
(|T9 ]i for q and p. By construction, due to the unitarity of U(t), the commutation relations for 
the system operators are preserved; also a direct verification is possible by showing that the 
quantum stochastic differential of [q{t),p{t)] vanishes due to d^Tt . Our choice of the field 
state will be such that the mean values of the noises Cq{t), Cp{t) are vanishing and this gives 
that the evolution equations for the mean values of q and p coming from (fT9l l are exactly the 
classical equations for an underdamped oscillator. This fact is a first justification of the choice 
(dill for the Hamiltonian and of the not usual connection (fThl l of position and momentum with 
the mode operator. 


2.1 The field state 

As field state we take the mixture of coherent states 

= E[|e(/T))(e(/T)|], /t(s) = 1(o,t)(s)/(s), (22) 

where / is a complex stochastic process with locally square integrable trajectories and E de¬ 
notes the expectation with respect to the probability law of the process /. In the argument 
of a coherent vector only square integrable functions are allowed, while the trajectories of 
the process / are only locally square integrable. So, we have introduced the cutoff T, rep¬ 
resenting a large time, which we will let tend to infinity in the final formulae describing the 
quantities of direct physical interest. As explained in 03 Sect. 3.2.1] this is a field analog of 
the regular P-representation for the case of discrete modes 0. In quantum optics, mixtures 
of coherent vectors with respect to true probabilities are interpreted as classical states. 

To represent the phonon bath 03 we take / to be a complex Gaussian stationary stochastic 
process with vanishing mean, E[/(f)] = 0, and correlation functions 

nfit)f{s)]=0, E[mf{s)]=:F{t-s). (23) 

Thanks to stationarity, the function F{t) is positive definite, so that according to Bochner’s 
theorem its Fourier transform 

/ -l-oo 

dt (24) 

-OO 

is a positive function, which we assume to be absolutely integrable, thus implying a finite 
power spectral density for the process. 

By this choice of the state we get that the noises (l20l i have vanishing means and symmet¬ 
rized quantum correlations given by 

^^{{Cq{t),Cqis)}) = ^^{{Cp{t),Cpis)}) 

= jra — [Sit-s) + 2ReF{t-s)], (25) 
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[d{t- s) + 2ReF{t- s)], (26) 


where ({C'i(i), Cj (s)}) := limr^+oo Trr {{Ci{t),Cj{s)} Let us stress that the noises 
appearing in system ( fT9] l cannot be arbitrary. They have to guarantee the preservation of the 
commutation relations [q{t),p{t)] = i by satisfying suitable commutations relations, equa¬ 
tions d2Tli in our case. Moreover, their symmetrized correlations must be compatible with 
their commutators, as it must be 


E 


52 


dt / dsh,it)h,{s) ^({Q(f),Q(s)} + [a(f),C',(s)]) > 0, 


VT > 0, for all choices of the test functions hi{t). Also this property is true in our case be¬ 
cause our noises and their correlations are an exact consequence of a unitary model and of the 
choice of a well defined state. This is not true in other proposals, where the positivity prop¬ 
erty above is not satisfied or divergences are introduced by not well defined approximation; 
see the discussion in lfT4l Sect. 3.3]. 


2.2 The reduced state of the mechanical oscillator 


Let po be the initial state of the oscillator. It is easy to see that the random reduced state 
Trr{C/(f)po® |e(/T))(e(/r) |C^(f)^} satisfies an usual quantum master equation with random 
coefficients. But this is not true for its mean, the reduced state 

p{t) = Trr{C/(f)po ® 0<t <T. 


By the properties of HP-equation there is no dependence on T as long as T > L In any case 
it is possible to characterize the equilibrium state of the system by solving the linear quantum 
Langevin equations (fT9l l and computing the first two moments of q{t) and p{t) for t -foo. 
The reduced equilibrium state 

Peq= lim lim Trr (p(0) OcTt^) (7(f)^} 

t->-l-oo T—s-l-oo 


turns out m to be a Gaussian state with {q)eq = {p)eq = 0 and 


(«")., = (p')., = — («.»+ 5 ). 

Wm V 


({9,p})eq = -^fiVeff+i), 
Wm \ 


where 


iVeff := 


7m r N{v) 

2tT (Win - vf 


diy. 


(27) 

(28) 


An important property of our model is that the energy equipartition in mean holds: 

('Z^)eq = (p^)eq- Moreover, the mean of the Hamiltonian (fTTI) turns out to be 


(77m) eq — 


7Veff+- 


3 Radiation pressure interaction 

We consider now the case of a mirror mounted on a vibrating structure and directly illumin¬ 
ated by a laser, so that it is subject to a radiation pressure force. One has to add a further inter¬ 
action term into the HP-equation suitable to produce a force proportional to the photon flux 
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in equation (I19bb for p. If we consider a well collimated laser beam and a perfect mirror, it is 
possible to represent the light by a single ray impinging on the mirror and reflected according 
to the laws of the geometrical optics. So, we take d — 2 and Bi{t) = H — H^, 

Ri = Om, S’!! = 1 as before; moreover, we add a further field i? 2 (f) = repres¬ 

enting the electromagnetic field, and we write Aem(f) = R 22 {t). A force proportional to the 
rate of photon arrivals means to have a term udAem(f) in (I19bl i. By comparing this expres¬ 
sion with ([Tol l with X = p, one sees that we need S 22 = S = and 512 = ^21 = 0, 

i ?2 = 0; 0 is a phase shift introduced by the mirror. So, the final HP-equation is 

dC/(f) = H^dt + (amdBl^{t) - aijdBth(t)) 

+ (5 - 1) dAe„,(f) - ^ a3„a„,dt|[/(f), (29) 


S = u e R, ()) e [0,27r), (7(0) = 1. 

From (doll one gets the relevant quantum Langevin equations 


dam(f) = 

(lWm+ ^'^a 4 n{t)dt ^/^dBth{t) +1TV^J dAem(f), 

(30) 

or, equivalently. 

dq{t) = HmP(()d( -1- dCq{t), 

(31a) 


dp{t) = - (n^nqit) + 7mP(f)) d( -b dC'p(f) -b t;dAem(f); 

(31b) 


The quantum noises Cq{t) and Cp{t) are given by (l20l l. The linearity of such equations allows 
for an explicit solution 


a^{t) = f e-('“”+^)(‘-^)dBth(s) 

+ lo (32) 

leading for the position and momentum Heisenberg operators to 

—7nii/2 f 4 . I 7m9 + , \ 

q(t) = e ' q cos ujmt H-sinwmf 

V 2a;in J 

+ f sinwm {t - s) dAem(s), (33) 

JQ 

U\ —7mt/2 f 4 2f2„i9 + 7mP . 

pit) = e ™ ' pCOSWnif-SinWmf 

V 2a;in J 

+ U J ^COSWm (f - s) - sinWm (f - s)^ dAem(s). (34) 
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3.1 Input-output relations 

We now consider the Heisenberg picture for the electromagnetic component of the field; 


BZ\t) = AZit) = U{t)^A,Zt)U{t). (35) 

By (fTSl) . (fTSl) we get the input-output relations 

dBZit) = 5(f)dSem(i) = (36) 

dA°"*(f) = S{t)^S{t)dA,Zt) = dAe„,(f). (37) 


Note that the number operator for the photons is not changed by the interaction with the 
mirror. 

By using (l3^ the scattering operator can be decomposed as the product 


Sit) = = ^o(t)Wth(^i; ll)Wem(0; 

where a system operator and two Weyl operators appear: 

«S'o(f) = exp I 


, 7ni7 “f , 

gcosWincH-^- sinwmr 


2Wn 


t—^ + oo 


)A’th(^t; 1) = exp I J £t{s)dBl^^{s) 


— h.c. 


>Vem(0; = exp i 


HmU 


2 ft 


e sinwm [t - s) dAem(s) f- , 


JO 


with 


iti>) = - 




iV{t)u){s) = Vis;t)u{s), Vu G L'^ 
Vis]t) = exp{iu2A(t- s)l(o,t)(s)} , /i(r) = 


n 


m _ Im. r • 

— e 2 sin ojrnT- 


(38) 

(39) 

(40) 

(41) 

(42) 

(43) 


The Weyl operator Wth(^t; 1) (l40l) is a displacement operator with function it (l42l i acting 
on the thermal component and >Vem(0; 14(f)) (HTt is a Weyl operator acting only on the 
electromagnetic component and characterized by the unitary operator V(t) (l43t . 


3.2 The field state 

Now the environment is described by a two-component field and its state must describe the 
phonon bath and the laser light. As field state we take the mixture of coherent states 

o'Jnv = E[|e(uT))(e(ur)|], ut(s) = 1(o,t)(s)u(s), m(s) = , (44) 

where / is the stochastic process described in Section IZTI and g describes a phase-diffusion 
model of a laser Qa, namely 

p(f) = woAO, Lp > 0; 
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W(t) is a standard Wiener process independent from the process /. It is easy to see that 


lim 

T^'+oo 


1 

7t 




g(t)dt 


2 . 


\M"Lp 

r 2 r, ? 

-^ + {g — OJq) 


so, the laser light has cari'ier frequency ujq and Lorentzian spectrum of width Lp. A possible 
generalization would be to take A —^ A(t), with X{t) a further stochastic process; this would 
allow to describe also amplitude fluctuations. 

With this choice of the state, for the thermal noises Cq and Cp we have vanishing means 
and symmetrized correlations (1251) . (l26l l. while for the electromagnetic held we get 


Tr {dSemW^Jnv} = Ac )‘df, Tr = |A|' dt, 

Tr{dB|^(s)dBem(i)o-Jnv} = 

Tr {dAem(s) dAem(t)o-J„v} = “ s) + |A|^ |A|^ dtds. 


(45) 

(46) 

(47) 


3.3 The equilibrium state of the mechanical oscillator 


Again, we can introduce the reduced state of the mechanical oscillator 

pit) = Trr{t/(f)po ® (^InvUit)'^}, 0<t <T, 

and the reduced equilibrium state 

Peq= lim liin Trr {(/(f) (p(0) 0 (7(<)1'} . 

t—>-+oo T—^+ca 

By working in the Heisenberg picture, from (l32l i- (l34t and the moments of the fields we get 
easily 


(p)eq — 0; (9)eq — 


4^ 


=: qo 


({g,p})eq = -^ Weff + - , 




(9^)eq -qlc = (4)eq = — ( Weff + ^ ) + 
Wn 


WV 

27m ’ 


(48) 


(49) 


where A^eff is given by (l28l l. By (l49l) the energy equipartition in mean holds again for the 
fluctuation part. Moreover, the mechanical mode occupancy is given by 

/„t \ ^raqlo .T I 44 |A|^^ 

\^m^m/eq rt _ — Aeff “t“ 


2lu„ 




we have also 


(®m)eq 






2Wn 


dWinl/m 





= ir 


4V 

4Wm 


Finally, it is possible to show that in the limiting case of constant phonon spectrum, 
N[v) Weff, and no phase diffusion, Lp ), 0, the reduced system state satisfies a Markovian 
master equation with Liouville operator 


/-[p] — t [T/miP] + 7m (Aeff + 1) | UniPOm n {®m®m)P} 


+ 7m/Veff ^amPOm “ ^ {omaln.p}^ + |A|^ (e‘”'^pe - p) . 
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The last term is new and describes the momentum kicks due to the scattering of photons. 
The other terms of the Liouville operator have the appearance of an usual generator for the 
dynamics of a mode in a thermal bath; however, the important point is that the link of the 
mode operator with position and momentum is not the usual one, but it is given by (flSl l. 
(HHl El Sect. 2.2]. 


4 Heterodyne detection 


To get information on the mechanical oscillator we can detect in various ways and analyse the 
light reflected by the vibrating mirror. In the balanced heterodyne detection scheme the light 
coming from our system is made to beat with a strong laser field (the local oscillator); the 
light impinging on the mirror and the local oscillator are produced by different laser sources; 
the stimulating laser frequency wq and the local oscillator frequency, say fj,, are in general 
different. Moreover, the phase difference cannot be maintained stable and this erases some 
interference terms. It can be shown lfT9l Sect. 3.5] that the balanced heterodyne detection 
scheme corresponds to the measurement in continuous time of the observables 


/(/r; t) 



dBem(s) 


+ h.c., 


(50) 


where a is a phase depending on the optical paths and > 0, represents the 

detector response function. In the Heisenberg description the observables become the “output 
current” 

=U(t)'‘I{^i;t)U{t). (51) 

By using (l36l l we obtain the explicit expression 


= J{t)+ h.c., J{t) = 


xe 


i(“+0) 


g-f (t-s)+i/iSgt 


(52) 


By the definition of /(/r; t) and the properties of U (f) (see the discussion at the end of Section 
[T]) we get [J(^; s)] = [/out(ft; t), s)] = 0, which says that the output current 

at time t and the current at time s are compatible observables. Note that to change /r means 
to change the frequency of the local oscillator, that is to change the measuring apparatus. In 
general /out (ft; t) and /out(ft^ s) do not commute, even for t = s. 

By the rules of quantum mechanics, once one has the commuting observables /out (ft; t), 
t > 0, and the system/field state po ® cr^^, the probability law of the stochastic process 
representing the output of the detection apparatus is obtained iniiisi. By taking the second 
moment of the output current the mean output power is obtained Gum and at large times it 
turns out to be proportional to 


P(p) = lim i [ (/out(ft; f)^)Tdf, (•)t := Tr {• po (g) crf„^} ; (53) 

T^.+oo 1 Jq 

the limit is in the sense of the distributions in p. As a function of p, P{p) is known as power 
spectrum. 

By using directly (I53I) . (fSTT i. (fSOl l. dJTl l. without computing the explicit expression of 
P{p), one gets easily the “total output power” 


1 

27r 


dp[P(p)-l] = 2\X\\ 


(54) 
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For sake of comparison it is interesting to have also the power spectrum of the input light; by 
setting K := + Lp, we get 


-Pin(/r) = lim i / (/(/i;t)^)rdt = 1+ ^ (55) 

T-S.+00 T Jo 

the final explicit expression in (l55T l is easily computed by using (l46l l and the CCRs. Moreover, 
we have immediately 

± f d^l[PM-l] = 2\X\\ (56) 

Jr 

Let us stress that the equality of the total input and output powers is essentially due to (iJTT i. 

4.1 Exact results 

The explicit expression of the power spectrum can be computed, as we shall show below. 
Firstly, (l5^ reduces to 


P(/r) = 1 + 2 lim {J{tyj{t))Tdt; 

T^+oo T Jo 


(57) 


then, we obtain 


^+oo 

P{p) = 1 + 4 |A|^ exp<{ 2 |A|^ Re J du “ 


{N, 


eff 


n 


m o 


r*+oo 


X Re / dfeO('^-‘"«)-t)*expi |A|" / ds 


^+oo 


dv'^h(t+s) _ 2 


+ / di/ 


+ 1 ) ' 


Wii 

Q-iv'^h(s) _ 

N{iy)e 


dTTOJm + {l^ - 


, (58) 


where h{t) is given in (021 and k = X + Lp. Let us stress that this is an exact result obtained 
from a unitary quantum evolution and the monitoring in continuous time of commuting ob¬ 
servables. Note that in this expression the thermal contributions, the terms containing N{y), 
and the electromagnetic contributions, the terms containing the function h, are completely 
interlaced. 


4.1.1 Proof of equations (l57l i and (l5^ 

Let us sketch now the proof of the previous formulae. By (l52]) we have 

(7out(ir; f)^)T = 2Re(J(f)^)T + (J(f)J(f)^)T + (J(f)^ J(f))T- 


By the presence of the limit in (152 . these terms contribute to P(/i) only with their large time 
behaviour. By using (l36l l. (I38l l- (l42 l. (l46l l we get 

^ ^^ 2 g 2 i(Q-i- 0 ) f ds [ dr 
Jo Jo 

X L(s;r)(Wth(4;l)>Vth(4;l))T(>Vem(0;L(s))>Ven.(0;L(r)))T, 
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Jo Jo 

X (>Vth(4;l)^Wth(4;ll))T(>Vem(0;l^(r))^>Vem(0;F(s)))T, (59) 


(J{t)J{t)^)T -1-2 k |A|" Re f ds T dr R(r; s)e—o)(s-r) 

Jo Jo 

X e-^(*-’')(Wth(4;ll)Wth(4;ll)t)(Wer„(0;R(s))Wer„(0;'l/(r))^)T. 

Then, one can check that limt^+oo limT-).+oo(-/(i)^)T = 0. Moreover, by using the com¬ 
position law ® for Weyl operators and (l42T i. (l43T l. we get 

Wem (0; R(s)) Wem (0; F(r)) ^ = Wem (0; R(r)) ^ Wem (0; Vis)) 
and, for s, r large and s > r, 

R(r;s)Wth(4;ll)Wth(^r;l)^ 

= v{r; s)Wth(4; a)^Wth(4; 1) exp { 2 ilm (f«|4)} ~ >Vth(^r; ll)^>Vth(4; 1). 

This gives (J(f) J(f)'l')T ~ 1 + (J(f)) 7 ' and, so, (l57l i is proved. 

Let us consider now (l59l i: recall that /i(rt) and V(s, t) are given in (l43T l and ft(s) in (l42li . 
Firstly, the electromagnetic contribution gives, for s > r, 

(>Vem(0; R(r))^yv’em(0; R(s)))t = expj ^ ^l/(M;r) R(m;s) - 1^ |A|^dM 

'^h{s—u) _ ^ ^—iv^h{r—u) _ ^ 


= a(5) a(r) exp< |A| / du 


~ |a(cx))| exp< |A| / du 


h{s—u) _ ^ 


^—iv^h{r—u) _ ^ 


with 




iv^h{u) 


- 1 


a(s) = exp-j 

Then, again for s > r and both large, the thermal contribution gives 
(Wth(^r;l)^>Vth(4;l))T 

¥sf + \Kf' 


= E 


exp i 2 ilm (/|4 - Ir) + ( 414 ) - 


= exp 


pr \ pr _ 

|fs('w)|^dM + J |4(w)|^duJ+ J ir{u) ^siu)du 


-3^. 


due*“"(4(u)-4(M)) 


~ exp< - 


i‘2NeS + 1 ) 


2uj„ 


2^^“ -^R 27ru;m ( 


^mJmN{i/)v^ cos iy{s — r) 


^ + (Wm — 

./» 47ra;„ 

By inserting these results into (l59l l and the expression found into (ISTT i. we get the final result 


(27Veff + 1) , f , 

= exp<--h / dv 


13 























4.2 Linear response 

When the laser light is used as a probe to get information on the dissipative oscillator, the 
beam can be taken to be weak, which means | A|^ small. In this case only the linear response 
is important and we can simplify (l58l) by considering only the “optical susceptibility” 


S(u) := lim 
|A|fO 


Pjf^) - 1 

|A|^ 


By (l58l) the weak probe limit gives immediately 


S(^) = dexp'l — 


(iVeff + i) r! 


m o 




X exp 


^+oo 


Re 






[{Njiy) + 1) e---* + iV(z/)e-" 
dTTWin + {iy - 


(60) 


Now, the power spectrum is 

P(M)=il + |A|'s(/r) (61) 

and we see that in this limit the thermal contribution is completely unaffected by the elec¬ 
tromagnetic one. So, we can use the optical probe as a mean to gain information on the 
mechanical occupancy spectrum N{y). 

To compute the time integral in (l60l l one needs to develop the last exponential in a power 
series. The result is much more clear when N{y) is slowly varying in a neighbourhood of 
Win of width 7 ni. In this case we can made the approximation N{v) ~ N{ujjn) in the last line 
of dhOl l: by (l28T l we have also By power expansion we get 


E(/r) ~ 4exp<^ - A—^i^ Re 


Wn 


(* + CO 




X exp-i - 


L L 


(A^(a;in) + 1) + iV(wi„)e-0‘"”+^)' 


(W^(Wni) + 2 ) 2 

= 2exp< —- — - V 

OJm 


00 m 


^ (A^(Win) + l)'^iV(Wni)’" •^(K + m7in) 

"" 6i±ipa)! + _ 2j) Wm]" ■ 

So, E(^) appears to be a series of peaks centred on ujq ± nwm and we write 


E(^) 2exp^ 


UJ„ 




(63) 




By reorganizing the sums we get the expressions of the various peaks and by integration their 
weights. 

• The peak centred in Wo: 


00 

no(M) = E 

j=o 


(iV(Wm) + l)'’iV(Wi„)'^ (n + 2 j 7 m) 


{jiy Adum 


4 


{^l - UoY 


(64) 
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here the term with j = 0 represents the elastic scattering of photons, while a term 
with j > 0 represents the scattering of a photon with exchange with the mechanical 
oscillator of j energy quanta Wm- The weight of the peak is 


^ I no(p)dp = ^ ■ 
i=o 


9_ / —^ , ,2 

1=0 uO 


1 (n^v^{NM + i)NM' 


4Wn 


(65) 


For N (wm) = 0 the previous formulae reduce to 
no(/t) = - 


X + ~ '^o) 

• The peaks centred in ujq — nuj^a, n = 1,2,..(Stokes scattering): 

\i+" 


tx / no(M)dAi = 1- (66) 

27r Jr 


n-n(M) = 


^2j-+n^4,-+2n(jv(^^) + + (2j + n) 7 n 


f=0 


j\{j + n)!22j+"a;S 


2i+n 


(K+(2i+n)7„) 2 

-2-^ + (p - Wo + nwm) 


(67) 

here the term with j = 0 represents the cession of a quantum Wm from the photon to 
the mechanical oscillator, while a term with j > 0 represents the same process plus the 
exchange of other j quanta. The weight is 


1 

27r 


3=0 


. / n_„(p)dp = (iv(wi„)+i)” V 

7b .-_n 

For N (wm) = 0 we get 

n_„(p) =- 


n‘^J+nv^3+2r^^NM + iyN{u3^y 


j\{j + n)!22J+"wS 


2j+n 


( 68 ) 


n (k + njm) 


n!2"w" 


(« + "7m) , ,2 

+ (/r - Wo + nwm) 


1 /■ n”?; 


n„,2n 


(69) 


(70) 


The peaks centred in wq + nWm, n = 1,2,..., (anti-Stokes scattering): 


il 


j!(j -f n)!22i+"wm^ 


2j-|-n 


(K-|-(2i-|-n)7m) . ,2 

^-7-^ -f (/r - Wo - nWm) 


here the term with j = 0 represents the cession of a quantum Wm from the mechanical 
oscillator to the photon, while a term with j > 0 represents the same process plus the 
exchange of other j quanta. Note that the weight turns out to be 


1 /'rr . f 

— / n„(/x)dAi = ———— 

27r Jr \ A^(wm) + 1 


1 

27r 


n_„(p)dp. 


(72) 


For iV(wm) = 0 we get n„(/r) = 0. 
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The asymmetry between Stokes and anti-Stokes scattering is the base for using the the 
optical probe as a device for thermometry at low temperatures. Indeed, we have 

4n_i(^)dM-4ni(p)dp ^ 

and this quantity can be estimated by the area under the curve of the experimental data when 
the peaks in ojq ± Wm are well separated from the elastic peak in loq, which means that the 
widths 7 ni and k = x + Lp are sufficiently small. The resolved-sideband thermometry is a 
technique already used in somewhat similar situations 101161. 
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